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A WEAK LAW OF LARGE NUMBERS FOR A LIMIT ORDER BOOK MODEL 
WITH FULLY STATE DEPENDENT ORDER DYNAMICS 

ULRICH HORST AND DORTE KREHER 


Abstract. This paper studies a limit order book (LOB) model, in which the order dynamics depend 
on both, the current best available prices and the current volume density functions. For the joint 
dynamics of the best bid price, the best ask price, and the standing volume densities on both sides 
of the LOB we derive a weak law of large numbers, which states that the LOB model converges to a 
continuous-time limit when the size of an individual order as well as the tick size tend to zero and the 
order arrival rate tends to infinity. In the scaling limit the two volume densities follow each a non-linear 
PDE coupled with two non-linear ODEs that describe the best bid and ask price. 


1. Introduction 

While limit order books have extensively been discussed in the economic and econometric literature 
for some years (cf. for example [3, 8| .12, ;21]), they have only recently gained increased attention by 
researchers in mathematical finance. One research objective is to specify a realistic discrete dynamics 
of a LOB which can be approximated by an analytically tractable continuous time model. This is 
achieved by introducing scaling parameters and passing to the high frequency limit, when the number 
of submitted orders gets large while the individual order size and the tick size tend to zero. Depending 
on the scaling assumptions the high frequency limit will either be deterministic as in a law of large 
numbers or be of (jump) diffusion type as in a functional central limit theorem. 

Deterministic high frequency limits for LOB models were derived by [15] and [11] . In [15] a weak 
law of large numbers is established for a limit order book model with Markovian dynamics depending 
on prices only. In nn the authors study a limit order book model, similar to ours but without any 
feedback effect, and derive a deterministic ODE limit using weak convergence in the space of positive 
measures on a compact interval. A diffusion limit for order book dynamics can be found in [6], where 
the top of the book is analyzed. The result was later generalized in [Tj. In m a high frequency limit 
for a one-sided limit order book model is derived under the assumption that on average investors place 
their limit orders above the current best ask price. The opposite case when orders are placed in the 
spread with higher probability is analyzed in [18] , where the authors use a coupling between a simple 
one-sided limit order book model and a branching random walk to characterize the diffusion limit, 
cf. also [22] . In the recent preprint m the limiting behaviour of an individual order position together 
with the best bid and best ask queue is studied and fluctuations around their fluid limits are derived. 

There is considerable empirical evidence (see, e.g. [31EUH] and references therein) that the state of 
the order book, especially order imbalance at the top of the book, has a noticeable impact on order 
dynamics. However, in the literature the order flow in most limit order book models either follows 
independent Poisson dynamics or depends on the price process only as in m- Exceptions to this are 

1991 Mathematics Subject Classification. 60F05, 90B22, 91B70. 

Key words and phrases, limit order book, market microstructure, high frequency limit, fixed point iteration. 

This research was partially supported by CRC 649: Economic Risk. We thank Moritz Greving for his assistance with 
the simulation results. 


1 



H]> wh ere Hawkes-type dynamics are used, as well as [E] and the very recent preprint []7|, in which 
the ergodicity of a general Markovian order book model is studied and the diffusivity of the rescaled 
price process in this general framework is derived. 

In this paper we adapt the model from [15] but use a different approach which allows us to deal with 
much more general, in fact fully state dependent Markovian order flow dynamics: the type of order 
(market order, limit order placement, cancellation), its size, and the price level at which the order 
is submitted can all depend on the current state of the limit order book, i.e. on prices as well as 
on standing volumes. This is different from [2], where the standing volume only influences the price 
dynamics, but there is no direct feedback to the order flow. As a result, unlike in m, the price process 
cannot be analyzed separately. Instead, we have to establish joint convergence of prices and volumes. 
The resulting scaling limit for a fully state dependent Markovian order book dynamics is the main 
result of this paper. 

Our main theorem states that when the number of submitted orders goes to infinity over a fixed time 
horizon, while the proportion of active orders, the tick size, and the individual order size tend to zero, 
the dynamics of the prices and the volume density functions converge to the unique solution of a non¬ 
linear coupled ODE/PDE system. To prove our main result we first construct a deterministic discrete 
non-linear approximation S^ to the random discrete order book dynamics S^ n \ This is done using 
a weak law of large numbers for triangular martingale difference arrays as in m, even though our 
method of approximation is different and more elegant, which allows us to handle this more general 
setting. In the next step we then construct an iteration towards the deterministic approximation for 
fixed n, denoted S^ ,m , and we prove that it approximates S almost uniformly. Afterwards it is 
shown that each iteration step in the prelimit converges as n goes to infinity to a continuous model S m 
solving a certain differential equation. Indeed, these models can be seen to be a fixed point iteration 
generated by a contraction mapping. The fixed point then gives a solution to our limiting coupled 
ODE/PDE system. 

For the ease of notation we have chosen to analyze only the buy side of the order book together 
with the bid price in most parts of this paper. However, if one defines the sell side and ask price in 
an analogous way, the result can easily be extended to a two-sided order book with order dynamics 
depending on the whole limit order book, i.e. on bid and ask prices as well as the order volumes of 
both sides of the book. The corresponding result for the two-sided LOB model can be found in the 
final section of the paper. Especially, making the distribution of order types depending on the bid-ask 
spread will ensure that the bid and ask price do not cross, cf. also [15] • Moreover, we assume that 
order arrival times are deterministic. However, one can easily generalize our main result allowing for 
randomly spaced arrival times by making use of the time change theorem as has been done in mm- 

The remainder of this paper is organized as follows: In Section [2] we define the dynamics of a sequence 
of one-sided discrete limit order book models, state our assumptions and the main result. We also give 
an example which satisfies all our assumptions. Section [3] is devoted to the analysis of the limiting 
coupled PDE/ODE system, while Section [4] contains the convergence proof of the discrete order book 
models to the high frequency limit. Finally, in Section [5] we state our main theorem for the two- 
sided limit order book model and conduct a simulation study, which shows how the state dependency 
can be used to equilibrate the buy and sell side volumes from an initially highly imbalanced volume 
distribution. 
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2. Setup and main result 


In this section we define for every n € N a model for the dynamics of a one-sided limit order book 
with tick size Ax^ n \ Later we consider the scaling limit of these models when the tick size and the 
impact of a single order tend to zero, while the number of order placements and order cancelations 
over a given time horizon [0, T] tends to infinity. Our modeling framework closely follows |l5j but we 
allow for a much more general dependence of order arrivals on the state of the book. Throughout, all 
random variables are defined on a common complete probability space (fi, T, P). 

2.1. The model. The dynamics of the order book in the n-th model is described by a cadlag stochastic 
process gW = (g( n )(t)) 0<t<T taking values in the Hilbert space 

- 2 / 


B:=MxL 


a 


|E ■— | a l| T 11CK2 11^2 ■ 

The state of the book changes due to arriving market and limit orders. In the n-th model there are 
[T/Au n 'J such events taking place at times 

4 ”^ := kAt^ n \ k = 1 ,... 4 T 


A fW 


where Adenotes a scaling parameter converging to zero as n —X oo and tg = 0. 

(n) 

The state of the book after k events is denoted S k ; . We put 

:= S<"> := for t £ [(<">, (&) , 

(n) (n) 

where B k and v k J denote the best bid price and the buy side volume density function relative to the 
best bid price ( relative volume density function ), respectively. To be precise, defining := j Ax^ 
for j G Z, n E N, 


L 


(n) 


S n ) 


(n) 

— 3—1 


(x)dx 


represents the liquidity available for buying at a price which is j E No ticks below the best bid price 
at time t), '. In order to model placements of limit orders inside the spread, the function v, . k E N, 
will be defined on the whole real line. We refer to the buy volumes standing at positive distance from 
the best bid price as the shadow book, cf. Figure [I] 


Quantity 


Shadow book 



-6 -4 -2 o +2 +4 +6 Distance from the 

best bid price in ticks 


Figure 1. The shadow book extends the relative volume density function to the right 


3 





















The idea of the shadow book is taken from m- The shadow book has to be understood as a tool 
to model the (conditional) distribution of the size of limit order placements inside the spread in such 
a way, that those placements extend the current volume density function of the visible book in a 
sufficiently “smooth” way to the right. The shadow book follows the same dynamics as the volumes 
of the visible book and becomes part of the visible book through price changes. The working of the 
shadow book and its interaction with the visible book will be further explained below; cf. Example 
E31 

At time t = 0 the state of the limit order book is deterministic for all n£N and denoted by 

4 n) = e«x L 2 (M). 

To state the convergence condition on the sequence of initial states we introduce for each n E N the 

(n) (n) 

translation operators T) and T_ , which act on functions / : M —> M in the following way: 

Ti n) (/)(•) :=/(• + AsW) , Ti n) (/)(-) :=/(•- A*W) . 

Note that the translation operator is isometric, i.e. for all / E L 2 , 


T[ n \f) 


L 2 


L 2 


Furthermore, let us fix some constant M > 0 throughout. 

Assumption 2.1. The initial volume function Uq"' 1 is a non-negative step-function on the grid {x^\ j E Z}, 
which is uniformly bounded by M and has compact support in [-M, M] for all n E N. Moreover, there 
exists a non-negative continuously differentiable function vq G L 2 such that 


(n) 

Vq -VQ 

(n) 


L 2 


= O 


i^Ax^ 


n) 


Also there exists Bq G M+ such that B q —> Bq. We denote sq := (Bq,vq) G E 


Remark 2.2. Note that Assumption 2.1 implies that vo also has compact support in [—M,M] and 
therefore, because vq G C 1 (M), its derivative must be bounded by some C > 0. Thus, 


T (n) _ j. 


X n ) 


L 2 


< 


< o 


->( n ) 


(n) 

Vq J -v 0 


L 2 


+ 


T 


(n) 


- 7 ) (vo) 


+ 


(n) 

VQ ~ v 0 ' 


(Ax^ n A + sup | Vq (x) | 
' ' a:SK 


Ax^ 


L 2 


L 2 

= O (Ax (n) 


L 2 


In our model there are three events that change the state of our limit order book. The (buy side) 
limit order book changes if: 

• (A): a market sell order of size equal to the current best bid queue arrives. In this case the 
best bid price decreases by one tick. Hence, the relative volume density function shifts one 
tick to the right. 

• (B): a buy limit order is placed inside the spread one tick above the current best bid price. 
In this case the best bid price increases by one tick and the relative volume density function 
shifts one tick to the left. 

• (C): a buy limit order placement of size at price level occurs. If < 0, this 

corresponds to a cancelation of volume. 

Here Av^ is a scaling parameter that determines the size of an individual placement / cancelation. 
We refer to market orders and limit buy order placements in the spread (Types A,B) as active orders. 
They lead to price changes. Cancelations and limit order placements (Type C) do not lead to price 
































changes. They are referred to as passive orders. The assumption that market orders match precisely 
against the volume standing at the top of the book and hence shift prices by exactly one tick is made 
for convenience and shows that our framework is flexible enough to allow for larger market orders. 
However, it is not unrealisitc: in an empirical study the authors of |10] found that in their data sample 
around 85% of the sell market orders which lead to price changes match exactly the size of the volume 
standing at the best bid price. The effect of a market order that does not lead to a price change is 
equivalent to a cancelation of standing volume. 


Event types are determined by a field of random variables ( d>[ n ' > ) taking values in the set 

V h J k,n £N 


{A,B,C}. The size and the price level at which an order placement resp. cancelation takes place 


are determined by a field of random variables 
tion. 




> Pk 


/c,nG No 


according to the following assump- 


Assumption 2.3. There exists a field of random variables ^ ir k 
interval [-M, M} almost surely and 


(n) 


fc,nS No 


taking values in the compact 


p■= 44 + for 7 r 


(n) 


g 


S n ) J n ) 

j-i’j 


Furthermore, there exists a field of random variables ( 
k, n G Nq. 


(n) 


such that 4 n) G [—M,M] for all 


fc,neNo 


The random variables 4"^ . k,n G No, determine the placement/cancelation price levels relative to 
the best bid price. If p j/ 1 = 4 > then the placement/cancelation takes place at the best bid price; 

(n) (n) 

if p k ; < B k , then it takes place deeper in the book; else it takes place in the shadow book. The 
shadow book interacts with the visible book through price changes which shift the relative volume 
density functions v^ n \ The following example illustrates the working of the shadow book. 

Example 2.4. Suppose that the kth event is a limit order placement one tick above the best bid price 
into the shadow book, i.e. 

4 n) = C , pt ] = B[ n) + Ax<"> and > 0. 

Further suppose that the (. k + l)sf event is a buy limit order placement in the spread, i.e. (p' k ^_ { = B. 
Then, 

4 % = 44 + Ax{n) = 4 n) + Ax(n) 

and for all x G [—Aa/ n ),0) corresponding to standing volumes at the current best bid price, 

4+20*0 = 4ti (* + A ^ (n) ) = 4 n) (*+ Aa;(n) ) + 

while for all x ^ [—Ax^ n \o), 

4+2 ( x ) = 4+i ( x + = 4 n) . 


Note that in general is (even conditionally) dependent on 44 one wan ts to avoid negatives 

volumes due to cancelations. The main contribution of this paper is that the conditional distribution 

of the random variables ( (jlfi ' 1 , , 7r[ n ' ) ) may depend on both, current prices and volumes. This 

V J fc.neNo 

extends m where prices are independent of volumes as well as [2j where only the distributions of 
price increments depend on volumes. 
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For every n G No and k = 0,1,..., [T/At^\ we define the u-field := a (^Sj U \ j < k'j. We 

will assume that for each n G N the state process S^ is a Markov process in its own filtration, 
cf. Assumption 2.5 below. 

To formulate the next assumption we need to introduce the space E' := {s = (B,v) G E : v G C 1 }. 

Assumption 2.5. 

(1) There are two Lipschitz continuous functions p A , p B : E —> [0,1] with Lipschitz constant L and 
a scaling parameter A p^ such that for all n G No and k < \T /A^ n )j, 


d n ) 


= I 


Ei n) ) = a 


SI 


(n) 


a.s. for I = A, B. 


(2) There are Lipschitz continuous functions f : E —>• L 2 , n G No, with common Lipschitz 
constant L > 0 such that for all k < [T/At^\, 




(n) 

k 


T' 

J i 


(n) 


a.s. 


and 


sup 

s€E 


/ ( ”'W(-) 


< M. 


Moreover, there exists a function f : E —>• L 2 such that 


sup 

s&E 


/ (n) M - m 


L 2 


= o 


(^Ax^ 


n) 


where /[s](-) : M —>• is continuously differentiable in x for all s G E’ with derivate 

being uniformly bounded in absolute value by M. 


Remark 2.6. Note that Assumption 2.5 implies that 


sup 

sGE 


_ j ' 


)(/ < " ) W)| ij = o(ax< 


n) 


The following example illustrates how our modeling framework allows for a dependence on the price 
dynamics and standing volumes. 


Example 2.7. Given any h G L 2 define 


ff<“> := 


v^ l \x)h(x)dx. 


(n) 

We may interpret Hi as an indicator for the volume standing at the top of or deeper into the book, 


depending on the choice of h. Set 


P 




_ , i ( o( n ) rr( n ) 

— 9 [B k ■ tl. 


I = A, B, 


for Lipschitz continuous functions g 1 : M 2 
exists L < oo such that 


[0,1]. Then by the Cauchy-Schwarz inequality there 


P 


s ; 


(n) 


P 


s. 


(n) 


< L 


g W _ g{n) 


Let us further assume that for all n G Nq, k < [T/ At^J, / = A, B, 


4 n) = i 


S \ n) , j < k) = A p (n, p 




s\ 


(n) 


a.s. 
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and that there exists a function m : M' 

E dx 


L 2 such that for all k, n as above 


Sj n \ j < k) = m 




(x)dx, 


where for each ( b , h) E M 2 the function m[b, h](-) is uniformly bounded with bounded support in [-M, M] 
and the mapping ( b,h,x ) i->- m[b,h\(x) is continuously differentiable with bounded derivatives. More- 


over, suppose that there is a field J of i.i.d. random variables with bounded density which has 
compact support in [— M,M] and set 

o A f or some e > 0. 

If the random variables 7r[ n \ and are conditionally on {S^ n \j < k} independent and ifoj^ 
is also independent of{S^ n \j < k}, then Assumption 2.5 is satisfied. 


2.2. Main result. We are now ready to define the full dynamics of the order book. For notational 
convenience we define for I E {A, B, C} and k, n E N, the event indicator variable 


:= lj 


(n) 

k 


and introduce the short-hand notation (/ = A,B ): 
pW- 1 [■] := A p (n V[■], pW’ B ~ A ; = p W’ B - p("). A , 

Definition 2.8. For each n E N the dynamics of the state process S ^ is given by 


An),A B-A B _ A ( n),B-A (n),B _ - (n),A 

P , P ■—P P , JLfc ^-k 


II 

'e 

J n ) 

s o 


£ 

(1) 



v i 

where 



A t(") 


i”> =-a+( t-a 1 - /) (»a) ifc’f +(h"> - /) (4A) ita + 

(n) 

<’(•) := ,,<;>)}(■)• 


To derive a law of large numbers we need to make the right assumptions on the scaling parameters. 
Our choice of scaling introduces two time scales, a fast one for limit order arrivals and cancelations 
and a comparably slow one for market order arrivals and limit order placements in the spread. 

Assumption 2.9. There exist constants co,ci,c 2 > 0 and (3 E (0,1) such that 

Ai^A p^ Au (n ) Aa;W 

lim - —f—. -= c 0 , lim . . = ci, lim - -* = c 2 . 

re->oo A t\ n ) n-A oo A A n ) rwoo (Af(™))^ 

TT./.o.g. we assume that Co = ci = c 2 = 1 m f/ie following. 

Remark 2.10. While it is very natural to assume that Av^ ~ At^ n ) for n —> oo in order to keep the 
total volume of orders in the limit order book of constant size, the assumption Ax^Ap^ ~ AfO) 
is not so standard. However, note that this constitutes the critical (and interesting) case. Indeed, as 
can be easily seen from the proof of our main theorem, assuming that Ax^Ap^ = o(Af^ n )) would 
lead to a constant price in the high frequency limit. Such a result can be found in m- 


The following weak law of large numbers is the main result of this paper. It states that the state 
process converges in probability to a deterministic limit that can be described as the solution of a 
system of non-linear differential equations subject to an initial boundary condition. 
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Theorem 2.11. Under Assumptions 2.1, 2.3, 2.5, and 2.9 there exists a deterministic process S : 
[0, T] —>• E such that for all e > 0, 


lim P 

n —»oo 


sup 

\0<t<T 


S (n \t)-S(t ) 



= 0. 


The function S = ( B , v) is the unique classical solution to the following coupled ODE/PDE initial 
boundary value problem: 

5(0) = s 0 , 

(2) dB(t) =p B ~ A [S(t)\dt, t G [0, T], 

v t (t,x) = p B ~ A [S(t)]v x (t,x) + f[S(t)](x), (t,x) G [0, T] x M. 


The rest of the paper is devoted to the proof of Theorem 2.11 In the next section, we first show with 
the help of a fixed point argument that the ODE/PDE system {[2]) does indeed have a unique solution. 
Section [4] contains the proof of convergence of the discrete models to the scaling limit. The difficulty in 
proving Theorem 2.1 1| comes from the non-local dependence of the coefficients on the whole function 
in (|2|. 

Example 2.12. Let us choose , n ,/' ) , and as in Example 


LO 




= u 


k > 11 k 
(n) N 


1 


M + l 


UJ 


(n) 


H n) (4 n> ) 


2.1 and set for some M > 0, 
AM . 


Then, choosing M large enough, the ODE/PDE system |i|) takes the special form 


5(0) = s 0 , 

dB{t) = p B ~ A [B(t), H(t)]dt, t G [0, T], 

vt(t,x) = p B ~ A lB(t),H(t)\v x (t,x) + fi[B(t),H(t)\(x) - f 2 [B(t),H(t)](x)v{t,x), (t,x) G [0,T] xR. 


Throughout the paper we will denote by C > 0 a generic constant that may vary from line to line and 
is independent of any index involved. 


3. The limit model 


In this section we prove existence and uniqueness of a solution to ([2]). First, we explicitly construct 
a solution via a fixed point iteration on a suitable Banach space to prove existence. Under the 
assumptions of Theorem 


2.11 


there exists a constant K < oo such that 


E S, 


(n) 


j( n ) 


< Bq 1 + T 


Ax( n ) A p( n ) 
AfM 


+ 




L 2 


+ T ■ 


Av^ 


f H [s] 


L 2 


< K 


for all k, n G N. We thus choose as our Banach space the space E of functions g : [0, T] —> E which 
satisfy 

sup \\g{t)\\ E <K 

te[o,T] 

equipped with the norm sup te r 0T i ||g(t)|| s . Uniqueness will be shown using a standard Gronwall 
argument. 
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3.1. Fixed point iteration in the scaling limit. To construct a solution to ([ 2 ]) we perform a fixed 
point iteration for the function F : E —> E defined through F : g t—G, where G = (G^G 2 ) and 
G 1 : [0, T] —»• M and G 2 : [0, T] —> L 2 are given by 

G'(t) = B 0 + [ p B ~ A [g(s)\ds, 

Jo 

G 2 (t, x) = v 0 (x + J P B ~ A b(s)] ds^j + J f [ 5 ( 5 )] ( K x + P B ~ A L 9(u )] ^ ds. 


Note that by dehnition and Assumptions 


2.1 


2.3 


and 2.5, G 2 (t, •) has support in the compact interval 


[-AT — T,M + T} for all t G [0, T], 

We define E' := {g £ E \ g : [0, T] -> £'}. 

Lemma 3.1. For fixed g € E' the function G = (G^G 2 ) satisfies 


G(0) — so, 
dG 1 (t) = p B ~ A [g{t)]dt , 

G t 2 (t,x) = p B ~ A [g(t)\G 2 x (t,x ) + /[fif(t)](x) V (f,x) e [0,T] x R. 

Moreover, in this case G G E' and there exist two constants J,L < 00 , which do not depend on g, 
such that 

\G 2 x (t,x)\<J, \G 2 (t, x)\ <L V (t,x) G [0, T] x R. 


Proof. It follows from the general theory of first-order PDEs that G solves the claimed PDE. Moreover, 
by Assumptions |2.1| and |2.5j 

|G 2 (t,x)| < sup | Vq (x) | + T- sup \{f[s])'(x)\ = : J < 00 , 

ckGM. #EM 

seE' 

\G 2 t {t,x)\ < \G 2 x (t,x)\ + \f[g(t)](x)\ <J + M =: L <00 


as well as 


sup ||G(t)|| E < B 0 + T + ||v 0 || L 2 +T • sup ||/[s]|| L 2 < K. 
te[o,T] sge 


□ 


tells us that it must indeed solve (PL 


If we find a fixed point of F which lies in E', then Lemma 3.1 
To do this, we will show that the function F is Lipschitz continuous to conclude with Banach’s fixed 
point theorem. In the following we write Fs,t(g) : = G 1 (t) and F V)t (g) := G 2 (t ) for G defined as above. 
Then: 

\F B ,t(g) ~ FB,t(g)\ < [ \p B - A [g(s)}-p B ~ A [g(s)]\ds <2L f \\g(s)-g(s)\\ E ds. 

Jo Jo 

Moreover, 

II- F v,t{g) - F v,t(g)\\ L 2 < 


VO ( • + / P B ~ A [<?(*)] ds )-vo[- + I p“- A \g{ S )] ds 

rt 


/ 0 


+ 


to 


f bO)] ( • + / P B A [&(«)] du)-f b(s)] ( • + / P B A [g(u )] du 


L 2 


„B-A 


ds. 


L 2 
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Since Vq is uniformly bounded with bounded support the mean value theorem along with Assumption 
2. 5| yields, 

v °(^ + j Q pB A [s( s )] ds ) ~ V ° (' + J 0 pB ~ A ds ) 2 - C J o ll#( s ) “ 9(s)\\ E ds. 

Similarly, as |(/[s]) / (-)| is uniformly bounded and has bounded support in [— M, M] for all s E E' and 
as / is Lipschitz continuous, 


/ b(«)] (■+ [ p b A b(“)] du) - f [s(s)] (■+ [ P B A [.g{u )] du 


L 2 


< C J \p B A [a{u)]-p B A [g{u)}\du+ \\f[g(s)]-f[g(s)]\\ L2 

< 2CL J s h( u ) ~9( u )\\E du + L \\9( s ) ~9 ( s )\\e- 

Hence, there exists some I\ > 0 such that for all t £ [0,T], 

\\Ft{g) ~ F t {g)\\ E < K f \\g(3)-g(8)\\ E ds. 

Jo 

Now, the space E is also a Banach space with respect to the equivalent weighted norm 

—at 


:= sup e~ at \\g (t)|| 


for any a > 0. Choosing a := 2 K we get 


o <t<T 


1 


and 


\F t (g) - F t (g)\\ E < K J e 2Ks \\g-g\\*ds < -e ZKt \\g-g\ 


\\F(g)-F(g)l<-\\g-g\ 


Therefore, by Banach’s fixed point theorem there exists a unique fixed point S of F. As noted 
above S solves Moreover, the sequence of continuous time models defined via S° = so and 

5 m+ i := F , m 6 No, converges to S. We have the following result. 

Theorem 3.2. The fixed point S solves the ODE/PDE system Q) and 


lim sup 

te [ 0)T ] 


S m (t ) - S(t) 


= 0. 


The following lemma shows that S m = yB rn ,v m j is Lipschitz continuous with respect to time. 

Lemma 3.3. There exists a constant L > 0 such that for all m G N and all s,t E [0, T\, 

S m (t)-S m {s ) < L\t — s|. 

E 

Proof. W.l.o.g. s < t. By the mean value theorem there exists some u £ (s,t) and for every x £ M a 
point u x £ (s, t) such that 


B m (t) - B m (s ) 


B m (u) 


{t-s) = 


B-A 


(t-s)<t- s , 


\v m (t, x) — v m (s, x)\ = \v™{u x , x)| fit — s) < L(t — s), 

10 


























where the last inequality follows from Lemma 3.1 Since the function v m (t, •) and its partial derivatives 
have compact support in [—M — T, M + T], 

\\v m (t) - v m (s )\\ L 2 = ||l[-M-T,M+T](-)|^ m+1 (^,a;)|(i-s)|| L2 < y/2(M + T) • L(t - s). 

Setting L := 1 + Ly / 2(M + T ) it follows that 


S m (t) - S m (s) <L\t-s\ V m G N. 


□ 


3.2. Uniqueness. In order to show uniqueness of the solution to ([2]) we assume to the contrary that 
there exists another solution S and define the shifted volume density processes 


v(t,x) := v ( t,x — J p B A S(s) 


ds 


v(t, x) := v ( t, x — J p B A [5(s)]ds 


Then 


rt 


V t (t,x) = v t [t,x- j p B A S(s) 

rt 

„B-A 


ds ) — p 


B-A 


S(t) 


v x [t,x- J p 


B-A 


Sis) 


ds 


= f 


S(t) 


x- / p 

Jo 


Sis) 


ds 


and similarly for v. Integrating with respect to t and using that u(0) = -0(0) = vq = u(0) = u(0) we 
lain from Assum 

v(t)-v(t )\\ L 2 < 


obtain from Assumption 2.5 and the mean value theorem, 

S{u) dii \ - f [5(a)] (y~ J P B ~ A du 


f S{s) 


p 


B-A 


ds 


L 2 


< 


< C 


f S(s ) - f[Sis )] 


+ c / 
u Jo 


P 


B-A 


S(u) -P B a [5(u)] 


du I ds 


Sis ) - Sis) 


ds. 


Moreover, as \v x [t, -)| is uniformly bounded by J with bounded support according to Lemma 
by the mean value theorem 

f f*t \ / ft 

\\vit, ■) - v it, ■)\\ L 2 = 


3.1 


again 


J^ P B ^ [Sia)] ds^j -v(t,.-J o P B A [5(a)] ds^j 


< ||t;(t,-) -vit,-)\\ L 2 + 

rt 


c 


c 


Sis)-Sis) 


Sis)-Sis) 


v(t,- - j^p 

rt 

ds + c' 


B-A 


rt 


P 


B-A 


1 0 


Sis) ds)—vyt,- — j p [S(s)]ds 

•t 


L 2 


L 2 


Sis) 


ds — 


[ P B ~ A [Sis)} 

Jo 


ds 


ds. 


Furthermore, 


Bit)-Bit) < [ p B ~ A [5(s)l -p B ~ A [Sis)} ds<2L f Sis)-Sis) 
Jo L J Jo 


ds. 


Therefore, 


Sit) - Sit) 


< c 


to 


Sis)-Sis) 


ds 


and the continuous version of Gronwall’s lemma, cf. Lemma A.3 implies that S = S. 

li 






















































































4. Convergence of the discrete order book models 


The goal of this section is to prove Theorem 2.11 Note that as opposed to |15| we cannot treat the 
price process independently of the volume densities because the conditional event probabilities and 
order placements resp. cancelations do depend on both, prices and volumes. 

In the following we set for all n G N and k G M, 

k 


(n) 


Especially, this means that (T 


->(»») 


T 


-i 


(/)(•) :=f(- + kAx^) 


= Ti n) . The following important Lemma deals with multiple 
applications of the translation operator. It will be used repeatedly in what follows. 

Lemma 4.1. If f € L 2 is a step function in the Ax^ -grid, i.e. f satisfies f(x) = f {lAx^^j for all 
x G \lAx^ n \ (l + l)A:z;( n )) and all l G N, then for all k, h G M, 


->(”) 


- T 


-i( n ) 


(/) 


L 2 


< (\k — h\ + 1) || (t^ — /) (/) 


L 2 


Proof. Because is an isometry, it is sufficient to prove the claim for h = 0. To do this we first 

consider the case k E N and claim that in this case even 


((r| n) ) fc - (/) ^ < k || (r| n) - /) (/) 


L 2 


Obviously, this is true for k = 1. Assuming that it is true for k — 1 we get 


1 ) o 


((rf)'- 7 ) (/) 

((rf))"" 1 -/) (/) 


((r| n) ) fc ' -i) (/) 
(A n) - /) (/) 


L 2 


< k 


( r W _ ^ (/) 
( r H _ ^ (/) 


L 2 


L 2 


and hence the inequality follows by induction for all k G N. Next for k G (0,1) either 

( T W) fc _ (/)(*) = (rj n) - /) (/)(*) or ((if))* - (/)(x) = 0, x G R+. 


Therefore, in this case 


Now take any k G M+. Then, 


(if))*-/) (/) ^<||(T| n) -/) (/) 


L 2 


(r| n) ) fc -i) (/) 


L 2 


((^)) LfcJ -/) (/) 

< Lfej ||( A ) -/)(/)|| i2 

< (fe+ 1) | (r| n) -/)(/) 


+ 


T 


M 


- T 


-.(«) 


L*J 


(/) 


L 2 


((4 n) ) fc W - J ) (/) 


L 2 


L 2 


Finally, the general case follows from the isometry property of the translation operator: 


(if))*-/) (/) 


L 2 


->(«) 


—/c 


(if))*-/) (/) 


L 2 


i- It 


i(«) 


-fc 


(/) 


L 2 


□ 
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4.1. A deterministic approximation of the discrete model for fixed n E N. Recall that by 
definition 


— 1 -n (n),B — A 

4"’ = • ,„ 0 


k -1 


*j n) J + Av^ ^ (T 

1=0 


— 1 -n (n),B—A 

(?T,) \ 7 


M, 


(n) 


For each n E N we are now going to define two approximations to the discrete model dynamics S^ n \ 
a deterministic non-linear approximation S A) in which the event indicator variables are replaced by 
their averages conditioned on the state of the approximating sequence and a random approximation 

—(tt.) 

S in which the event indicator variables are replaced by their averages conditioned on the random 
state of the original state sequence. 

More precisely, we define for each n E N the process S A) through 

k -1 


B^ ] := B^ n) + A x ( n) Yp {n) ’ B ~ A 


l=o 


S 


(n) 


and 


v (n) _ ^ T |n)^ J=0 [ i J („(")) + A „(“) Y (Tj n) ) Ei=j ^ J (/("> 

1=0 


s. 


(n) 


for k = 0,1, 


A t(") 


~o( n ) 


Moreover, we define a second process S through 


k -1 


Rr := Rj n) + A x^Yp {n) ’ B ~ A 


(n) 


1=0 


and 


^ , Y^ fc_ 1 — A [c( ?1 )l , . k—1 „k — l (n),B — A [c( n )l , 

_(n) ;= ^ r |n )^ J=0 [ 3 \ ^(nP j + ^n) ^ M ( /(n) 

1=0 


5 


(n) 


for k = 0,1, 


Ad n ) 


In a first step we are now going to show that the sequence S (n> approximates S^ n \ The proof uses 
a weak law of large numbers for triangular martingale difference arrays, which can be found in the 
appendix. Subsequently we show the desired convergence of S ( n ) to the discrete model dynamics S^ n \ 

In what follows C > 0 denotes a generic constant that may vary from line to line and is independent 
of any index h, i, j , k,l,n. 


Theorem 4.2. Under the assumptions of Theorem 2.1 1\ for all e > 0, 


lim P sup 

n-s-oo l o <t<T 


S( n \t) - s {n \t) 


> e =0. 


Proof. For the bid price component we have 

k -1 


4 n) - ^ 


= Ax (n) 


v 1 W’- 8 _ |( n )A _ Jn),B 


1=0 


p- 


S. 


(n) 


+ P 


(n),A 


A n ) 
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By definition the random variables 


y-.M := /\ X H f i[ n )’ B _ jWd _ p (n),B 


: (n) 


+ P 


(n),A 


S 


(n) 


, 3 < 


A t( n ) 


n e N, 


form a triangular martingale difference array with respect to ). If we can show that there exists 

a > \ such that 


sup 

j<T/At 00 
neN 


/ E 

y. (n) 

2 


j 



(AiW) 


2a 


< oo, 


then Theorem A.l will imply that 


sup 


k< 


- A t( n ) 


k -1 

i=o 


(n) 


= o(l) in probability. 


Indeed, this follows immediately from Assumption 2.9 with a := ^4^, because 


E 


(n) 


< 4 (Ax^') -E (p {n)A 


(n) 


+ P 


( n ),B 


s. 


(n) 


) < 4 (Ax( n) ) 2 Ap^ < AC (a 


i+P 


Next we consider the volume component: 


(n) _(n) 

vi -Vi 


L 2 


< 


1 ti (n),B—A v^fc-1 (n),B —A ["oMl \ , 

j — (Ti n )| i= 0 P ^ ] (vq^ 


L 2 


+ 


Au (n) ^ (V_ 

i=o 

+ At,("> £ 

i=o 


k —1 — 1 -n (n) ,B — A 

(n)^ ^i=i i 


jVf( n ) __ jM 


5 


(n) 


L 2 


E fc — 1 -n — A 

*=3 * 


— I T 


. Y^fe-l ..(nJ.S-A re( n )] \ 

-,(n)^i=j p J ] f f (n) g(n) 


L 2 


Let us first deal with the second term. Due to the norm invariance of the translation operator this 
term equals 


k -1 




„ VV 1 II v 
(n)\ 4^i=o -i 


3=0 


— 1 it ( n),B-A 


\ ^ 

~sf 

) 



/ 


L 2 


The variables 


X^ n) : = 


•sr^j — l -n (n),B — A 

:= Au (n) ‘ (A/j n) - /W 


(n) 


form a triangular martingale difference array with respect to As 


M 


(«) 


is bounded by 


M/ Aan”) according to Assumption 


2.3 


E 


X 


(n) 


L 2 


< 



— 1 ii A 


«’) 

to 

II 

to 

m 



L 2 



L 1 


•E 


Mj 


(n) 


2 Au (n) ) 

< AM 2 - -— r — . 

Ll ~ Ax( n > 
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Therefore, for a := 44 > \ we have 


sup 

l<T7At( n ) 

neN 


E 

x (n) 

2 


3 

L 2 


(AfW) 


25 


< oo 


and hence by Theorem A.l 


sup 

fc 4ziWj 


k -1 

w 

1=0 


= o(l) in probability. 


L 2 


Regarding the third term note that by Lemma 4.1, Remark 2.6, and Assumption 2.9 
fc-i 


Au (n) ^2 


3=0 

k-1 


n (n)Ai=j H 


- T 


, —1 „(n),S —A [o( n )l ' 

^ I 6 * J 


(/ (n) 


5, 


(n) 


3 


< 


3=0 

k-1 

Au (n) ^ 


n ( n )< B ~A (n),B—A I c( n )l 

(n)Ai=j L Pi J 


- / j (/ (n) 


L 2 


1=0 




< C sup 
j<k -1 


fc—1 

^(n),B-A _ ^(n),B-A 

i=j 
k-1 

yy _ p (n),B-A s (n) 

i=j 


(n) 


+ 1 


+ 1 Az (n) . 


I. 2 


(rj n) - /) (/ (n) sj n) ) 


L 2 


Thus, we may conclude as above for the price component that the term converges to zero in probability 
uniformly in k < T / At^ n \ The convergence of the first term in the above decomposition follows 
analogously. □ 


Theorem 4.3. Under the assumptions of Theorem 2.1 1\ for all £ > 0, 


lim P sup 

n=>-oo l 0 <t<T 


S^\t)-S^ n \t) 


>£ \ = 0 . 


Proof. We have 

s (n) (4 n) ) - s (n) (4 n) ) 


< 


S (n) _ gin) ^ ||^ + I gin) ^ _ 5 ( n ) ^ 


tin) 

l k 


According to Theorem 4.2 the first term converges uniformly to zero in probability. In the following we 
use the Lipschitz continuity of the p 1 s and f^s formulated in Assumption [2A] to derive an appropriate 
upper bound for the second term, which allows us to apply the discrete Gronwall lemma in order to 
prove the assertion. We start again with the bid price component: 


B ( k ] -B < n) 


= Ax (ri) 


k— 1 

yyn),B-A 

1=0 


(n) 


_ p{ri),B—A 


in) 


k-1 


< Ax^ Ap( n hLj2\\Sj n) ~ S\ n) 

l=o 


Moreover, using Lemma 4.1 and Assumption |2.5| 

k—l —1 ( n ),S-Ara( n )l _ 

-,(n)\^i=j p J (An) q(n) 


A 

l=o 


^j( n ) gp) — f ^ S (n) ^ 


fc-i 


< Av( n >L^||^ n) ~s) n) 

L 2 1=0 


15 




















































































Similarly to the proof of Theorem 4.2 using Lemma 4.1 and Remark 2.6 we can derive the inequality 

„w,b-aL(»)] 


^rji(ri)\^i=3 p r* J _ frj,{n)y^i=i P 


f(n) 


S 


(n) 


L 2 


< c 


fc-1 

£ 

1=3 


P 


(n),B—A 


J(n) 


P 


(■ n),B-A 


S 


(n) 


+1 Ax (n) 


fc—1 


< C Ax (n) Ap (ri) 2L ^ Sf’ - 5> n) + CAx [n) . 

11 E 

i =0 


2.1 


Relying on Assumption 


computations to get a Sinn 


instead of Assumption 


2.5 


we may replace 


ar estimate for the initial volume term. 


5 


(n) 


by Vq L ' > in the above 


Finally, putting everything together we have by Assumption 2.9 

< 


~( n ) =( n ) 


L 2 


. \^k — l (n) ,B — A [ c( n )l sr^k — 1 (n),B — A [ c( n )l \ 

T (n)\Zj=0 R [ S J \ _ / T |n)\E,=o R L 5 ,' J \ / (n) 


L 2 


+ At,<"> £ 

j=0 

fc -1 

+ Av( n ) 

3=0 


Ti 


TT^k — 1 (n),B — A [c( n )l 

(n)\Uj J 


^y( n ) s^ s {n) ^ 


L 2 


. x=k — 1 (n),B — A [ci 71 )] , . V^fc — 1 (n),B— A [oMl \ „ 

Wp i=j P 4 J _ / r (n)^5J i=j R J f f(n) 


s. 


(n) 


L 2 


fc -1 


fc -1 


fc-1 


< C Af (n) ^ p} n) - s[ n) + Ax (n) + A f (n) ^ At (n) ||sf n) - s[ n) + Ax (n) 

\ z=o i=o L 2=0 

fc-i 

< CAt (n) ^ |s£ n) - s[ n) + CAx (n) . 

z=o s 

Therefore, for some sequence (a n ) converging to zero in probability and for all k < A ^ n) we get the 


following uniform estimate by means of the discrete Gronwall Lemma A.2 

o( n ) ~Q{ n ) 

O, — D i. 


g( n ) _ g(n) 


+ 


oi 72 ) o( n ) 
'-’fc ~ D fc 


fc -1 


1=0 


<a n + CAt^ p, (n) - s l n) 

< a n + a n CAt ("> Y ^- l+1 CM{n) < Ca n . 


fc-i 


1=0 


□ 


4.2. Almost uniform iteration to the discrete approximation. In this section we approximate 
for each n E N the model S^ iteratively by an if-valued sequence of limit order book 

models. To this end, we define for each n £ N a function FC) ; E -A- E via 


F (n) : g i—^ G [n) = (g^\g^ with G {n) (t) ■- ^ (n) s 


= Gj, if te 


(n) An) \ 

l i ’ T j+i) 


where G^ := and for k G N, 


u B,fc •“ ^B.fc-l + P 


,B—A 


(n) 

4-1 
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as well as 


» . 

),k 




k -1 


+ Au (n) Y 

3=0 





with 


■■= > ( f f) ■ 


We will write F^ n \g) := G'[. n) for G'j^ defined as above in what follows. 


Using analogous arguments as in the proof of Theorem 4.3 one can find a constant K > 0 such that 

~ 4+1® || E s AA " (n) E lift - ftlU + AA * ,n) f ® on k < » e N. 

3=0 


Now as in Section 3.1 we define a weighted norm on E via 


II0II** : = su p e 


—3Kt 


0 <t<T 


E > 


which allows us to get an estimate for the weighted norm with Lipschitz constant less than 1 up to an 
error of order As for all n E N, 

k 


( 0 ) — F y ' l> 
r k+l\y) r k +1 


(n) 


< 


KAvW Y e3Kt ^ } h ~ 9\L + KAxW 


= KAv {n) ■ 


3=0 

e 3K(fc+l)At( n > _ 


,3KAtW 


- 1 


-■\\g-g\\.* + KAxW, 


there exists by Assumption 2.9 an Nq 6 N such that for all n > Nq 

(9) - A i?,l 


—3 Kt^ 

e ^i k+1 


KAv^ 

< - -l|fl- 


, mu — oil + KAx^ < — llo — g|L + KAx^ 

E ~ e 3KAt( n ) _ i lla y|1 ** — 2 y a 


and therefore indeed 






V n > Nq. 


W.l.o.g. we take JVo = 1 in the following. For each n,m £ No we define a new discrete time model 
g( n )> m v i a S^ ,0 (t) = 4 ” } as well as 


c(n),m -\-1 


^?( n ) 



Theorem 4.4. 


lim sup 5 (n) ’ m (t)-5 (?l) (t) 

m ~ >co o <t<T 
n=fO o u — — 


: 0 . 


Proof. First, note that by definition is a fixed point of F^ n \ i.e. 

p(n),m p(n) Qm pH ^(*»)^ = 5 W V m <£ N. 


m times 
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Hence, making use of the above computations we deduce that for every m, n G N, 


Cj(n),m _ cj(n) _ p(n),m ^g(n),0^ _ p(n),m 


< 


< 


< 


v 7 fc=0 v 7 


sup 

2 1 0<t<T 

m— 1 


M_SH( t ) +2KAxW 


/i \ "i- 1 _ 

f-J A' + 2/fAi W . 

Now the result follows from the equivalence of the weighted norm and the norm sup 0<t<T ||-|| E . 


□ 


Especially, Theorem 4.4 implies that 


lim lim sup 

m—^oo n —yoo q <t<T 


sto’ m (t) - sW(t) 


= lim lim sup 

E n—y oo m —¥oo q<£<T 


- s( n \t) 


= o. 


4.3. Convergence of the discrete iteration to the continuous iteration. The goal of this 
section is to prove the following result. 


Theorem 4.5. For all m£tf, 


lim sup 

n-s-oc 0 <t<T 


S^’ m (t) - S m (t ) 


= 0. 


Proof. To prove Theorem 4.5 we proceed by induction. Obviously, for m = 0 the claim holds by 


Assumption 2.1 Now assume that the claim holds for m and consider the (m + l)-th iteration: First, 


we show the convergence of the bid price process. Writing the integral as a limit of Riemann sums we 


have by Assumptions 2.1 and 2.9 


[t/AiMJ 


B m+l (t) = lim Hi" + lim Ax^Ap^"'p 

n .—yon u n—yon * 


(n) B ~ A 


S m ( y 


3=0 


As p B A and S m are both Lipschitz continuous by Assumption 


2.5 


and Lemma 


3.3 


this convergence 


is uniform in t E [0, T]. Moreover by Assumption 2.5 


sup 
o <t<T 


Ax^Ap^ (p B ~ A 

3=0 


S m It 


i(n) 


z ^ AxWA p™ B—A 

5 r - Atw ' 

J ~ At( n ) 

Ax (n) A 

< 2LT •- . -• sup 

AR n ) 0<t<T 


S m ( ^ 


P 


„B—A 


„B—A 


(n),m 


P 


s. 


(n),m 


5 m (t) - 


which converges towards zero as n —> oo by Assumption 2.9 and the induction hypothesis. 

We now show the convergence of the buy side volume density function step by step. W.l.o.g. we only 
prove the convergence of the order placement / cancelation term. The convergence of the term involv¬ 
ing the initial volume density function follows by analogous arguments. First note that a pointwise 
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Riemann sum approximation gives 

ft r -I / pt 

x+ p 


f S m (s ) 


/o 


B-A 


S m (u) 


du ) ds 


|i/AtWj 

lim At {n) V / 

n—> oo ^^ 

J=0 


, t (n) 




B—A 


S m (u ) 


du . 


To show that the convergence also holds in L 2 observe that for all t < t G [0,T] with |i — t| < At^ n \ 
making use of Assumption |2.5j Lemma 3.3, and the mean value theorem, 


/ 5 m (t) 


+ / P 


B-A 


S m (u) du) — f S m (t) 


+ / P 


B-A 


S m (u) 


du 


L 2 


< 


/ 

"s m (i)' 

-/ 

(i)' 

+c r 

L 2 Jt 

P B ~ A 

>»] 


5 m (u)l du < LL\t - t\ + C\t - t\ < CAt^. 


Therefore, the above convergence does indeed hold in L 2 , uniformly in t G [0, T] . 
Second, by similar arguments 


Lt/At (n) J 

Ad") 

3=0 


(Tj n) ) 


d^ A ‘ WJ p ( " ) ' B “ A [s'i n),m ] 


< sup 

j<t/At( n ) 


sup 

zeR 


(/ [s m (*$ n) )]) , w 



'’)] (■ 

+ f p b ~ a l 


du'j 

/ |_t/At (n) J 


ft 


\ 

Ax (n) JZ P {n) ’ B ~ A 

g(n),m 

- / p b ~ a 

An) 

S m (u) 

du 

\ i=j 


Jt\ ’ 

J 


/ 


L 2 


L 2 


< C sup 
j<t/ At( n ) 


Aaj^ApW ^ p - 


B—A 


s 


(n),r 


L- 


B — A 


S m (u) 


du 


and as for the price component this term converges to zero uniformly in t G [0, T], 
Third, note that by Lemma |4.1| and Assumption 2.5 we have 


[t/At(">j 

f T H\A =j p \ /j 




A J2 

> (4 n) y 

-/ (n) 

g(n),m 


3=0 






< T sup 

j<t/At( n ~> 


S m t v 


X n ) 


-f (n) 


S m ( 


L 2 


+ 


/W 


S m t 


X n ) 


-f 


(n) 


L 2 

S. 


(n),m 


L 2 > 


< o(l) + L sup 

j<T/ At(") 


- s^ ,m (t^) 


which converges to zero by the induction hypothesis. 

Therefore, we proved that uniformly in t G [0, T] the term 

Lt/Ah»J Y^L*/ A *^ n b n (n),B-A\a(.n),ml 

A t<”) 53 ^ I 

3=0 

converges in L 2 towards 

r 


(f 


(n) 


(n),m 


+ i, p 


B—A 


S m (u ) 


du ) ds. 


□ 


19 


































































































4.4. Proof of Theorem 2.11| . Finally, let us put the partial convergence results proven in the 
previous subsections together to prove the convergence of the discrete limit order book models S 

an Mi = Mi(e) such that for all m > Mi, 


to S. For this fix e > 0. Then there exists by Theorem 

S m -S 


3.2 


e 

E < 6’ 


Also, by Theorem 4.4 there exist M 2 = M 2 (s) and N\ = iVi(e) such that for all m > M 2 and n > N \, 

^(n) _ g(n),m 


£ 

E < 6’ 


We set Mq = Mq(s) := M\{e) V M 2 (e). Then for all n > N±, 


S M ° - S 


+ 


g(n) _ a(n),Mo 


£ 

E < 3 


Furthermore, Theorem 4.5 yields the existence of an N 2 = N 2 (Mo,e) = N 2 (Mo(e), e) = A^e) such 
that for all n > N 2 , 

g{n),M 0 _ gM 0 


£ 

E < 6’ 


Hence, for all n > Nq = A r o(e) := Ni(s) V Wfs), 


S (n) _ g 


< 


g(n) _ g(n),M 0 


+ 


g(n),Mo _ gM 0 


+ 


s M ° - s 


£ 

E < 2' 


Finally, by Theorem 4.3 there exists for every 5 € (0,1) an IV 3 = N^(£,5) such that for all n > N 3 , 

gW _ g(n) 


< x ) > 1 -S. 

E 2/ 


Setting N = N(e, 5) := Nq(e) V A^e, 5) we conclude that for all n > N, 


( g(«) - g^< £)>!-£ 4 ^ 5 (n) - S > e) < 6 . 


As £, 6 > 0 were arbitrary, this completes the convergence proof. 


□ 


5. The two-sided limit order book model 


For the ease of notation we have concentrated on the one-sided LOB model in the previous sections. 
It is however straight forward to generalize the proof of Theorem 2.11 from a one-sided to a two- 
sided LOB model, in which both sides are modeled in a similar manner. In what follows, we briefly 
formulate the assumptions and dynamics underlying the two-sided model and state the corresponding 
convergence result. The proof will be omitted as it is works in exactly the same way as the proof of the 
one-sided model discussed in detail in the previous sections. To illustrate the usefulness of the general 
state dependency, Subsection |5.2| contains some simulation results of the two-sided LOB model. 


5.1. Setup and convergence result. In this section the state of the order book in the n-th model 
is represented by 

S (n Ht ) := , t € [0,71, 

where A^ n \t) describes the best ask price, B^ n \t) describes the best bid price, v^\t, •) represents the 
sell side volume density function and v^ n \t, •) the buy side volume density function (both in relative 
coordinates) at time t. Here, the volumes of the real order book correspond to v^ n \t,x), x < 0, and 
n \t,x), x > 0 , while the extension of v^ n \t) to the positive halfline and the extension of Va 2 \t) to 
the negative halfline will be understood as the respective shadow book of each side and are used to 
model the distribution of order placements inside the spread. 
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The event variables will now take their values in the set { A , B, C, D , E, E 1 } with 

• A = market sell order 

• B = buy limit order inside the spread 

• C = placement or cancelation of a buy limit order 

• D = market buy order 

• E = sell limit order inside the spread 

• F = placement or cancelation of a sell limit order. 


Here the buy side events A,B,C are defined in the same way as before, while the sell side events 
D , E, F are symmetric copies of the respective buy side events. Hence, the dynamics of the two sided 
order book can formally be written down as 


= Si n) := ( 


A 


(n) 


H 
^ a,k > 


v: u, B 


(n) 
k > 


( n ) 
%k 


t E 


An) An) \ 

L k > L k +1 J > 


where for k = 1 ,..., [T/At^l , 


4 ° = 

B\ n) = 


A H i + 1 (n),D- E Axin ^ 

B+ lWf- A AxW, 


(n) (n) 

= 




) (-'ll,) + iff 


_j 


I i|(n),F ( n ) ^ 

+1 fc-r 4») 6 

3 


•i". 'i)}^ 


.(") T (") 
o "H+ij 


(■), 


(n) 

V b,k = 


4ti + Aif (h”> - /) (<$_,) + h-’r 4 


j-’( n ) _ j 


A®W Wfc -i 4^ 


,(n) 

a.k— 1 


,(n) 

6,fc-l 


Besides the scaling assumption 2.9 we make the following assumption on the two-sided limit order 
book model in analogy to Assumptions 2 . 1 , 2.3, and 2.5 before. 


Assumption 5.1. 

(1) There exist A$,Bq E M+ such that A^ -A- Aq and B^ -A- Bq. Moreover, the initial volume 
density functions v^j and vj^ are non-negative step-functions on the grid {x^\ j E Z} ; which 
are uniformly bounded by M and have compact support in [-M, M] for all n E N. There exist 
non-negative continuously differentiable functions v a ,o,Vb,o E L 2 such that 


( n) 

v lo - v h 0 


L 2 



i = a,b. 


(2) The random variables 



k,n£ No 


and 



k,n£ No 


take their values in the compact interval 


[-M, M] almost surely. 

(3) There are Lipschitz continuous functions p A ,p B ,p D ,p E : E x E —> [0,1] with Lipschitz constant 
L and a scaling parameter A p^ such that for all n E No and k < \T/At^\, 


k {n) 


= I 


s 


(n) 


j < k) = A p 


Hp 1 


S-, 


(n) 
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a.s. for I = A, B, D, E. 

















(4) There are Lipschitz continuous functions : E x E —>• L 2 , n G No, with 

Lipschitz constant L > 0 such that for all k < [T/Aifo'J and I = C,F, 


common 


h n) 

and 


S , 


(n) 


H " AxW E ( w * } S 

' jiEZ l l y j 


(n) 

k 


sup 

sGE 


sf, j<k I o.». 


/{”’[*](•) 


< M. 


Moreover, there exist two functions fc , f F 'ExE^-L 2 such that for I = C, F, 


sup 

sGE 


/i n) H-//W 


L 2 


= o 


(W n) ) , 


where each //[«](■) : M —>• [— M,M] is continuously differentiable in x for all s G E' x E' with 
derivate being uniformly bounded in absolute value by M. 


In the following the norm on E x E is defined as ||(si, S 2 ) WexE := II s i|Ie + for all si,S 2 G E. 

Moreover, we set sq := (Ao, v a ,o, A), Vb,o) G E x E. 


Theorem 5.2. Under Assumptions 5.1 and 2.9 there exists a deterministic process S : [0, T] — > E x E 
such that for all s > 0, 


lim P sup 

n->°o \ 0 <t<T 


S( n \t)-S(t) 


ExE 


> £ = 0 . 


Moreover, S = (A,v a , B,v b ) is the unique classical solution to the following coupled ODE/PDE initial 
boundary value problem: 


5(0) = s 0 , 

dA(t) = p D ~ E [S(t)\dt, t G [0, T ], 
dB(t) = p B ~ A [S{t)]dt, t G [0, T], 
d t v a (t,x) = p D ~ E [S(t)\d x v a (t,x) + f F [S(t)](x), (t,x) G [0, T] x M, 

d t v b {t,x) = p B ~ A [S(t)\d x v b (t,x) + f c [S(t)]{x), (t,x) G [0,T] xR. 


In order to guarantee that the bid and ask price, the spread, and the volume density functions are 
non-negative, the following assumption has to be satisfied. 


Assumption 5.3. 

P B (s) = 0, P D (s ) = 
p A (s) = 

p (d n) < = c ) = 


0 V s G {(A, v a , B, v b ) G E x E : A = Bj , 

0 V s G {(A, v a ,B,v b ) eExE: B = 0} , 

0 V n G N 0 , k = 0,1,..., |T/Ai (n )j - 1, 

0 V n G N 0 , k = 0,1,..., [T/At ( foj - 1. 


Lemma 5.4. If in addition Assumption 5.3 is satisfied in the statement of Theorem \5.2j then both 
volume density functions, both prices, and the spread between the best ask price and the best bid price 
are non-negative. 
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5.2. Simulations. This subsection contains a simulation study of the two-sided limit order book 
model introduced in the previous subsection. It provides an example of the usefulness of the general 
dependence structure that is allowed in our model, but was not covered in the previous work [ITS] . 

For this let us fix h > 0. We define the imbalance factor in the n-th model at time t k ' via 

Im M , = VolBidP 

k ' VolAsk[ n) + VolBid[ n) 


with 


VolAsk jf } := / v%l(x)dx, 
Jo 


VolBid^ : = 


' —h 


v l n k ( x ) dx ■ 


(n) 

Moreover, in the n-th model the spread at time t, is defined as 

Sp^ :=A^-B^. 


There is strong empirical evidence, cf. e.g. [0| and [23], that the probability whether the next price 
change is upwards or downwards depends on the imbalance of the bid and ask queues at the top of the 
book: if the imbalance is high, i.e. the standing volume at the top of the buy side is significantly higher 
than the standing volumes at the top of the sell side, then prices are more likely to move upwards. 
And conversely, an imbalance factor which is close to zero will increase the probability of the price 
moving downwards. Thus, we choose 


4"> = A 
4”* = D 

An) 


= B 
= E 


Sf\ j < kj = A p (n) (l - ■ exp , 

Sj n \ j < kj = A p( n hm^ ■ exp , 

Sj n \ j < k^j = A p (n )/ m 0d • (l - exp , 

Sj n \ j < k'j = Ap (n) (l - Im^'j ■ (l - exp 


Moreover, for simplicity we suppose that placements outside the spread are always of size 50 (this 
assumption can easily be relaxed) and are located relatively to the best ask resp. bid price according 
to a Gaussian distribution, i.e. 

1 


4 rt) = = 50, 4 n) G dy 

4 n) = F, 4° = 50, 4 n) € dy 


Sf,j<k) = 


s { ;\j<k ) = 


v / 32'/r 2 

1 

v / 32'/r 2 


1 - A p (n) ^j (l - Im^ exp ( -y 2 ) dy, 
1 — A p( n ^ Im ^ • exp (— y 2 ) dy. 


We note that if integrated over the whole real axis the above two terms sum up to (l — A/2. 
Furthermore, cancelations are supposed to be proportional to the current volume. Hence, we choose 
for all x < 0 , 


= C, G dx, 7G dy 


S ) 


(n) 


j < k) = 


4A 


vi ,’{y) 


■ [-*>(»),»] (l) ^b (' _ ApW ) Im A ' exp dxiy 


IjH = p, ^(n) ^ ^(n) g dy 


S ) 


(n) 


j < k ) = 

1 


( n ha,\ ^ \/ 327T 2 


a,k 


(: v ) 


(x) 


^1 — A ^1 — Im^^j exp (—y 2 ) dxdy. 


If we integrate the above two terms in both variables and sum them, we get again (l — A/2. 
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For the simulation we choose the following parameter values: 


n = 250, a = 0.8, h = 0.55, T = 8. 

We start with a limit order book that has a severe imbalance at time t = 0: standing volumes at the 
ask side are many times higher than standing volumes at the bid side. However, our specifications 
about the conditional distributions of incoming orders made above guarantee that the imbalance at 
the top of the book will disappear after some time. While Figure [2] shows the evolution of the whole 
visible limit order book plotted at different points in time, Figure [3] only shows the evolution of the 
bid and ask price together with the cumulated volumes at the top of the book. 



Figure 2. The evolution of the limit order book volumes 


The evolution of prices is influenced by the spread as well as the imbalance factor. As can be seen from 
the first graph in Figure [3| first the spread is the dominating factor (due to the exponential) and forces 
the ask price to decrease faster than the bid price. However, after time t ~ 3 both prices decrease 
simultaneously keeping the spread almost constant. This downward shift of the mid price is caused 
by the imbalance of standing volumes at the top of the book and therefore mimicks the findings in 
the literature very well, cf. e.g. [3]. Also the second graph in Figure [3] shows that cumulated volumes 
at the top of both sides of the limit order book converge. In this particular simulation study the buy 
side volume approaches the sell side volume because we have chosen the size of order placements much 
greater than the size of average cancelations for the initial volume density functions. If placements 
were supposed to be of a much smaller size, for example 10 instead of 50, the opposite effect could be 
observed, i.e. the sell side volumes at the top of the book would decrease to approach the buy side 
volumes at the top of the book. Of course, much more general random and even state dependent 
choices of the order sizes are possible and would lead to even more interesting dynamics. 


Last but not least, Figure [2] shows that the discrete order book dynamics can indeed be well approxi¬ 


mated by a smooth function and is hence supportive of Theorem 5.2 
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Figure 3. Evolution of prices and of volumes at the top of the book 


Appendix A. 


A.l. A weak law of large numbers for triangular martingale difference arrays. The following 
weak law of large numbers for triangular martingale difference arrays relies on moment estimates from 
|20| and is shown in M- 

Theorem A.l. Let (yV , k = 1,..., n; n E N) be a triangular martingale difference array taking values 
in a real separable Hilbert space such that 

sup (n 2 "E|y£| 2 ) < oo 

k<n 
ne N 

for some a > ^. Then for all e > 0, 

lim P sup 

n—>oo \ m<n 



A.2. Gronwall lemmas. For reference we cite the discrete and continuous version of Gronwall’s 
lemma which can for example be found in [0]: 


Lemma A. 2. Let (y m ) m >o, ( fm)m>o, and (g m )m> o be nonnegative sequences. If 

m—1 

Vm< fm+^2 9kUk V m, 
k =0 

then 

m—1 

— 1 

fk9k^= k+19j ■ 

k =0 


Lemma A. 3. Let a, /3, u be real-valued functions defined on some interval I = [a, b], a < b, a, b, E M. 
Assume that jd and u are continuous and that the negative part of a is integrable on every compact 
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subinterval of I. If (3 is non-negative and if u satisfies the integral inequality 

rt 


then 


u(t) < a(t) + f /3(s)u(s)ds V t e I, 

J a 

u(t) < aft) + J a(s)/3(s) exp f}(r)dr^j ds, t€.I. 
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